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I. INTRODUCTION 

In quantum field theory, time and space play quite dif- 
ferent roles. One consequence of this is the production of 
particles in the presence of moving boundaries, which is 
particularly evident in the predicted response of the vac- 
uum electromagnetic field to a dielectric medium in mo- 
tion: nonuniformly accelerated mirrors should radiate pQ, 
and a process analogous to that of Hawking emission [2 
has been predicted to occur when a localized change in 
the refractive index moves through a medium at a speed 
exceeding the phase velocity of light in the medium [3]- 
5 . Indeed, it is surprising is that there should be effects 
of this type that occur for dielectrics in uniform motion. 
Nevertheless, one such effect is that of quantum friction, 
whereby two dielectric plates in relative lateral motion 
are predicted to come to rest due a production of pairs 
of surface excitations out of the vacuum state [6]. 

The physics of these phenomena been the subject of 
recent controversy, with some debate concerning the re- 
ality of quantum friction [TrfTT]. as well as difficulty in 
quantitatively interpreting the experimental observation 
of the aforementioned Hawking like process Q21 H3] . It 
seems like part of this discord arises from a lack of agree- 
ment about the theory that should be used to describe 
quantum electromagnetism interacting with a moving di- 
electric. Aspects of the quantization of electromagnetism 
in moving media have been considered before, initially by 
Jauch and Watson [2], but as far as the author is aware, 
there has been no canonical approach that includes the 
effects of dispersion and dissipation [33] . Here we de- 
velop a canonical theory that includes these effects and 
we apply it to some simple physical scenarios. 

The point of view taken here is that dispersion and 
dissipation are not merely complications that somehow 
superimpose onto a 'clean' case where a dielectric 
can be represented by a refractive index, n that is 
independent of frequency, but are in fact fundamental 
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to the interaction between macroscopic objects and 
the quantized electromagnetic field. After all, the 
fluctuation-dissipation theorem tells us that for linear 
dielectrics we can always understand the vacuum elec- 
tromagnetic field as being produced by currents within 
the dielectric that are in proportion to the square root 
of the dissipative response times h |15j . To develop a 
canonical theory of electromagnetism that is consistent 
with this perspective, we use a method similar to that 
of Huttner and Barnett [16] . where a fictional reservoir 
is added to the electromagnetic Hamiltonian to account 
for both the absorbed field energy and the dispersive 
response. Philbin has recently presented a modifica- 
tion and generalization of this procedure that applies 
to media that can be described by any e(x, u) and 
/i(x,ui) that satisfy the Kramers-Kronig relations [17j . 
Philbin's approach was extended within [18j, where the 
Lagrangian necessary to describe moving media was 
derived, the quantization of which we now investigate. 



II. CLASSICAL ELECTROMAGNETISM 
INTERACTING WITH A MOVING MEDIUM 

To begin, we give the form of the action necessary to 
derive the classical equations for the electromagnetic field 
interacting with a moving medium, as presented in [18] . 
This is given as the integral of a Lagrangian density over 
space-time, S[A^, X u , Y u ] = J d 4 xJzf, where 

The part of the Lagrangian associated with the degrees 
of freedom within the electromagnetic field is of the usual 
Lorentz-invariant form 

^ F = ^[E 2 -c*B 2 ] (2) 

where the fields are written in terms of the potentials 
(ip, A) as E = -Vip - A and B = V X A. The elec- 
tromagnetic field is coupled to two continua of oscilla- 
tors, & Y u that represent the collective degrees of 
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freedom of the medium. The coupling of the field and 
the oscillators is chosen so as to reproduce the experi- 
mentally measured linear susceptibilities of the dielectric 
from the equations of motion. For a dielectric moving 
with velocity V, the part of the Lagrangian describing 
the oscillators is modified from that of [17] due to the 
Lorentz transformation of the time coordinate 



jsf - 1 

L ./0 



duj < 7 



X^ + iV- V)X 

+ (X u , -r Y 



2 v 2 



(3) 



where 7 = (1 — V 2 /c 2 ) -1 / 2 . Additionally, the transfor- 
mation of the field amplitudes means that the coupling 
between matter and electromagnetism involves tensors, 
a EB & ck BE , in addition to those of [17] 



doj \a EE ■ X^ + a EB ■ Y wj 



/>oo 

+ B • duj [cx BB ■ Y u + a BE ■ X u ] . (4) 
Jo 

where the a tensors are typically functions of u>, x, and t, 
even if they are only functions of frequency and position 
in the rest frame. 

Considering the specific case where the medium is 
described in terms of a scalar permittivity and per- 
meability in the rest frame, and the motion is V = 
Vx, the coupling tensors take the form; a EE (x,ui,t) = 
Aa(x\uj); a BB (x,uj,t) = A/?(cc', w); a EB (x,uj,t) — 
j/3(x' ,w)l3 X V/c 2 ; and a BE (x,uj,t) = ~ja(x',u)t 3 X 
V, where a(x',uj) = ^2ujlm[x EE {x' , uj)]/tt, f3(x',u}) = 
y/2ujlm[x BB (x' , u})]/tt, and A = diag(l,7,7). The sus- 
ceptibilities, Xee & Xbb are those of the medium in the 
rest frame, and x' is the rest frame spatial coordinate: 
x' = "/(x — Vt); y' = y; and z' = z. As is well known 
(e.g. [H]), this modification of the coupling tensors turns 
an isotropic medium into an anisotropic one when it is in 
motion. Furthermore, if the medium is inhomogeneous 
in the rest frame, then it will be time dependent in the 
laboratory frame. Only when the motion is along an axis 
where the medium has translational symmetry will this 
not be the case. 
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B (6) 



with the Y u obeying the same equation after the substi- 
tutions E «-» B and X -f-> Y. Throughout the paper we 
shall most often only explicitly write out the quantities 
associated with the X u and the polarization, P, but al- 
ways with the understanding that the remaining results 
can be obtained using this substitution. 

The equations of motion for the electromagnetic field 
are found to be the usual macroscopic equations 



V • (e E + P) 
1 dE 



V X B 



c 2 dt 



= 
= Mo 
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~dt 



V x M 



(7) 



where P — duo [a EE ■ X u + a EB ■ Y u ], and M takes a 
similar form, but with the aforementioned substitutions. 
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FIG. 1: The electromagnetic field (E,B) is coupled to two 
continua of oscillators (X U ,Y u). The strength of the cou- 
pling is given in terms of the set of ex. tensors within Q, and 
the interaction of the field and these oscillators is such that 
the correct constitutive relations emerge from the equations 
of motion. In our case we are investigating a dielectric in mo- 
tion along x, that is allowed to be inhomogeneous only as a 
function of x« . 



A. Solutions to the equations of motion 

The solutions to the classical equations of motion of 
|TJ) form the basis of the quantization procedure, and we 
will expand the field operators in terms of these functions 
in the next section. 

The equations of motion for the continua of oscillators 
can be found from the usual formulae, 



d 
dxf* 



d^X, 



dXZ 



(5) 



With the simplification of translational symmetry 
along the direction of motion, x (see figure [TJ), the so- 
lution of ^ becomes 

^ ( , Q] _ cxUx l{ ,u)-E + o^ E (x l{ ,u;)-B 
" l ' 11 ' > ~ (u-fl-- it})(u + irf) 

+ X Qu ,(k,x [l ,il) (8) 

the tilde indicating a Fourier transform with respect to 
x, and t. In (|8| we have introduced the following new 
quantities: X 0uj = 27r[<5(f2_ — u)hx„ (k, ten) + <5(0_ + 
io)h* x (— sen)], O-t = 7(fi ± V • k), the vector xn is 
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the coordinate in the y — z plane, and r\ is an infinites- 
imal quantity that specifies the retarded solution of the 
driven oscillator equation (|6). The expression given by 
^ should be understood in the limit 77 — > and be evalu- 
ated in terms of poles and principal parts. The (classical) 
amplitudes, h Xuj may be freely specified, and represent 
the motion of the oscillators that is not driven by the 
electromagnetic field. 

Inserting ^ into we identify the electromagnetic 
susceptibilities as follows 



Xbb(^-) 



a EE (uj) • aLH 

'0 



a EB (ui) • ajLH 



■00 



(uj — f2_ — irj){ui + il_ + irf) 
a EE (uj) • a^ E ((j) + a 
(lu — 0_ — in)(u) - 



ir))(uj + 



■duj 



*eb(, 



(uj — 0_ — irf) [ui + Sl_ + ITj) 



(9) 



Notice that due to the additional terms depending on V 
within ([3]), the Doppler shifted frequency, £!_ appears 
within the susceptibilities. Furthermore, once Q has 
been expanded in terms of poles and principal parts, it is 
found that the real and imaginary parts of these suscepti- 
bilities are related by the Kramers-Kronig relations [TH] • 
Using ^ and the counterpart expression for the Y u , 
the & can be eliminated from The Maxwell 
equations then become those of macroscopic electromag- 
netism in a moving dielectric, defined in terms of the 
susceptibilities given in ([9]), and with 'free' current and 
charge densities proportional to h Xuj & h Yu . The for- 
mal solution of these equations may be obtained through 
finding an appropriate Green function. We thereby find 
the following expression for the electric field 



E(k, xu , fl) = ifl / d 2 x\ l G(k,x h x\,,n) ■ j (k,x\,,SV) 



where the Green function, G(k, xi\, xi, f2) is the solution 



of (B8|. The quantity appearing as a current in (10 1 is 



given by 

j (k, ac|| , fi) = -iflP (k, X|| , fi) + V x M (k, x\\ , f2) 

where we have defined, Pq — Jg £&^[q?ee * Xq^ ~t~ QIeb ' 
Y 0u] ], and V = kx + Vn. To complete the specification 
of the dynamics of the system, we can use the Maxwell 
equation, V X E = iilB, and (10) to determine the 
magnetic field. 



B. Classical Hamiltonian 

As the Lagrangian given in (JT) is local and contains 
only first order time derivatives, there therefore exists an 
associated Hamiltonian which we now derive. 

For the electromagnetic field the vector potential is 
the only dynamical variable and the associated canonical 



-e Q E 



momentum density is given by 
OA 

while the oscillators have the associated momenta 



n 



x„ = 



dX u 



= 1 I [X U1 + {V- V)X U ]. (11) 



Due to the fact that electromagnetism has a gauge sym- 
metry, the scalar potential has a canonical momentum 
equal to zero, 11^ = dJz?/d(p = 0, and is not a dynamical 
field. We must therefore interpret the first of ^ — which 
is not a dynamical equation — as a constraint that deter- 
mines ip after some choice of gauge has been made. 

From these quantities, we can form the Hamiltonian 
density 



J4? = U A ■ A 



n 



x^ 



which we give as the sum of two parts 

the first of these is 
1 



2e 

and the second is 



{u A + py - (V X A) -M 



(12) 
(13) 



(V x A) 

(14) 



n 



x^ 



-^Xt 



-2V-(V § 

+ (X L 



•x u yu Xu 



did 
(15) 



(10) Expression (14) was obtained through neglecting a term 



equal to a divergence (that will integrate to zero in the 
Hamiltonian) and then setting <^V • Ua = 0. This is 
the aforementioned constraint on ip: V • [eoE + P] = 0. 
Finally the Hamiltonian, H = H F + H R is obtained from 
the integral of ( 13 ) over all space. Notice that when V 



is set to zero we obtain the Hamiltonian of [171 . 



The effect of the motion of the medium in (14 15) is 
twofold. Firstly within (14) the canonical field momen- 
tum, 11^ is coupled to the Y u by the velocity, and A 
is similarly coupled to X^. This is due to the Lorentz 
transformation of the field amplitudes, and means that 
the medium appears as a magnetoelectric (i.e. the elec- 
tric polarization responds to the magnetic field as well as 
the electric field). 

Secondly, and more relevant to this discussion, the os- 



cillator amplitudes within ( 15 ) are coupled to the asso- 
ciated canonical momenta by the motion, via the term 
AH = - f °° V • (V <g) X u ) ■ n x „dw + (X u -> Y u ). 
This contribution has the physical meaning of the mov- 
ing medium responding to the Doppler shifted frequen- 
cies of the electromagnetic field. Indeed, if we want the 
susceptibilities ([9| to contain f2_ rather than fl then we 
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have no choice but to include AH in the expression for 
the energy of the system [18]. Yet this term is peculiar, 
and causes the Hamiltonian to lack a lower bound: the 
spatial dependence of X u can be made arbitrarily sharp 
so as to reduce AH to an arbitrarily negative value with 
no change in the rest of the Hamiltonian. Therefore, in 
quantum theory we should expect that the zero particle 
state will not be the lowest energy state of the system. 
Indeed, when we come to describe the system with quan- 
tum mechanics, we will find that even at T = OK, it 
is possible to extract energy from the vacuum fluctua- 
tions of the electromagnetic field outside of a uniformly 
moving dielectric. In appendix [C] it is demonstrated that 
AH is an accounting device that arises from enforcing a 
non-zero and uniform velocity that disappears once we 
include the centre of mass of the dielectric as a dynamical 
variable. Throughout this work we take the simplest case 
where we suppose that the velocity of the dielectric has 
been fixed to some constant value by an external force. 

As a final comment, we note that in reality the energy 
of a moving dielectric interacting with the electromag- 
netic field cannot be reduced to — oo. Spatial dispersion 
would become relevant at some large enough magnitude 
of AH, and provide a lower bound. Mathematically this 
would be evident through a dependence on V (X) X u 
in the rest of the Hamiltonian. This is an important 
approximation that does not seem to be mentioned 
within the existing literature. However, although this 
warrants further investigation, here it does not alter our 
conclusions so long as we are careful, it just means that 



in reality for a large enough velocity the Hamiltonian 
can be reduced to a very large negative value rather 
than — oo. We can avoid this complication for the time 
being so long as we ultimately avoid expressions that 
depend on wavelengths of the electromagnetic field that 
are comparable to, or smaller than the atomic spacing. 



C. Hamiltonian along a classical trajectory 

As in the vacuum theory of quantum electrodynam- 
ics, it will be advantageous to expand our field operators, 
(A(x, t),TlA(x, t),X u (x, t), IIx^ {x, t), . . . ) in terms of a 
certain basis. This expansion will be carried out in terms 
of the classical solutions given by ([§ and ([10]). Explicitly, 
the classical amplitudes hx u & hy^ that act as sources 
for the electromagnetic field (e.g. see ( |Xo| )) will become 
operators, analogous to the a(k) and d T (fe) of the vac- 
uum theory |20j . As we will be expanding the space-time 
dependence of the field operators in terms of these clas- 
sical solutions, we can use this to simplify the form of 
the Hamiltonian. This simplification is equivalent to the 
usual procedure of Fano diagonalization used in [161 117) . 

It turns out that the most straightforward route to 
simplify the Hamiltonian is to begin with the part as- 
sociated with the continuum of oscillators, H R . After 



applying (11), this is 



H R = / (fx,. 



dk f dQ 



; -i(n+n')* 



du}(oj 2 — f2_|_0' + 



x u (k, n) ■ x u {-k, of) + Y u (k, n) • Y u (-k, n r ) 



(16) 



Into (16) we insert (|8| and its magnetic counterpart, to 
which we apply the results, 2£l 2 _ — — = 



III. QUANTIZED HAMILTONIAN 



nipulations this leads to 



270(Q_ -O,), (Al ), and (B8). After some lengthy ma- 



/r f°° 2 
d 2 x\\ I — / 2uiuj+\hx, (aj|i,fc)| duj 
J 2-7T Jq 7 

+ (X -> Y) - H F (17) 



where we have introduced lo± = 7(0; ± Vk). The total 
Hamiltonian is therefore equal to an integration over 
the absolute squares of the amplitudes, hx~ w & hy 1U] 
multiplied by which can be either positive or 

negative. As in the previous section, it is clear that for 
certain choices of the amplitudes, it is possible for this 
classical Hamiltonian to take an arbitrarily negative 
value. 



To quantize the Hamiltonian given by the sum, (17) 



H F , we take the usual approach of quantum field theory 
and let the expansion coefficients become operators. In 
this case the expansion coefficients are given by the hx^ 
and hy^ , and we perform the following substitution 



C E (x\i, fc,7w) 



^ M (aJ|i,A;) -> 



hy^{x h k)^\l— C B (X||,fc,7w) (18) 



where it is assumed that the C operators are bosonic 

C x {k,x h Q),c{,{k',x[,Q') = 

27rl 3 <5 AA /5(Q - tt')5(k - k')5^\x\\ - x\ { ) (19) 
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The (normal ordered) Hamiltonian is then 



H = 



E 

A — E,B 



XII 



dk 
2^ 



C A ( 7W ) 

where for brevity we have suppressed the arguments, k 
and x» within the C\ operators. As expected, the lack 
of a lower bound for the classical Hamiltonian carries 



over to the quantum case, the eigenstates of ( 20 1 having 
eigenvalues from — oo to oo. 

It is shown in appendix [B] that this specification of 
operators and commutation relations is consistent with 
the canonical commutation relations that must hold 
between the fields and their conjugate momenta. It 
might be objected that we have substituted ( 18 ) in 



the final classical result (17), which itself assumes that 
hx^ & h* x commute. However, the Hamiltonian is 
quadratic in all variables, so in doing this we have only 
neglected a constant term which is equivalent to an 
unobservable phase factor in the overall wave function, 
and can therefore be dropped. 



IV. 



DETECTOR COUPLED TO THE 
ELECTROMAGNETIC FIELD 



The Hamiltonian ( 20 ) is interesting, but at present it 



is simply that of [T7j, described within a different frame 
of reference. To demonstrate the physical effects of a rel- 
atively moving dielectric, we must couple another system 
to the electromagnetic field. As a model we consider a 
point -like detector at a fixed position, x , with an in- 
ternal variable, x, and some internal energy states which 
are the eigenfunctions of a Hamiltonian, Hp. 

The detector (atom) is coupled to the electromagnetic 
field in the usual form dictated by minimal coupling [3T] . 
The interaction terms are then transformed into a gauge 
invariant form through a unitary transformation of the 
operators, O — UO'W where in the dipole approxima- 
tion, U ~ exp(f x ■ A(x a )) [21 US]. This then gives 
the following familiar expression for the full (interaction 
picture) Hamiltonian |34j 



H = H D -d- E(x Q ,t) + H m +f 



(21) 



where d — ex is the dipole operator that acts on the 
internal state of the atom, Hm+ f is given by ( 20 1 , and 



the electric field operator, E(x ,t), is expanded as in 
(B2|, with the coefficients (B6). 



We consider the Hamiltonian associated with the inter- 
nal states of the atom to be in the form of equally spaced 
energy levels, Hjj = hua^a, and the interaction to be of 
the form, Hj(t) = -in ■ E(x , t)(ae" iw * - a 1 
the a & are the usual raising and lowering operators, 
and k is a constant vector. Note that in the presence of 
a lossy dielectric, the atom is coupled to the full electric 
field, rather than just the transverse part pMl |2"5] . 



It is not possible to introduce a set of orthonormal 
number states to use with the C operators. In vacuum 
QED one would introduce periodic boundary conditions, 
and the dispersion relation would then restrict both the 
wave-vector and the frequency to discrete values so that 
a set of such states could be introduced. In our case 
the dissipation within the medium removes the dispersion 
relation between the frequency and the wave-vector for 
excitations within the medium, and this option is not 
open to us. A straightforward alternative is to use a set 
of orthogonal states with infinite norm, defined over a 
continuum, as in [26] . However, here we use the approach 
of |27j . and work in terms of a new set of operators for 
which number states can be introduced, labelled by four 
integers (l,m,n,p), 



dk 
2^ 



d 2 xn 



dn 



x 4>i, m Ai,p{k,x\\Vi)Cx^{k,x\\,Vt) 



where we have expanded the normal mode operators in 
terms of unit vectors, e a : C\ = J2a e <rC\,<j, and chosen 
the orthonormal basis functions, 4>i, m ,n.p 



— I d*xu I dilt 
2n 



'h,m 1 ,n 1 ,p 1 ( Pl 2 ^ m 2,n2,: 



&l\l2 &m\rri2 ^n\ri2 ^piP2 



that form a complete set, 



<^m,n,p(fc, £C|| , ^)<t>l, m , n ,p( k '^ X \\ > °') 

l,m,n,p 

= 2TrS(k - k')5 (2) {x\\ - x\)5{tt - VL 1 ) (22) 

These new operators then satisfy the familiar commuta- 
tion relations 

[c\ ucn (h,mi,ni,pi), cx 2 ,a 2 {h, m 2 ,n 2 ,P2)} 

— ^AiA2 3<Ti&2 ^1^2 ^m\vn,2 ^n\ri2 ^p\P2 

and we can use this fact to define number states, 

1 



\ n \,a : i,j,k,l) — 



c[ji,j,k,l) |0) (23) 



If we suppose that the detector is initially in the ground 
state and the combined system of field and dielectric 
is in the zero particle state, |0) |0), then to first or- 
der in |k| the only possible transition for the detector 
is into the state |1), and the field-matter system into a 
state, |lA,<T,»,j,k,l)> We therefore write the complete wave- 
function as, 



), where \ip) = |0> <g> |0) + ^ C\,a,ij I *,i(*)|l> ® |1a w ,m> 

The usual results of time-dependent perturbation the- 
ory can now be applied [28 . The rate of change of the 



G 



quantity C\,a,i,j,k,l(t) is, 



■^(lA,«r,w,fc,i|®<l|if7(t)|0)®|0) (24) 



Inserting the expression for the interaction Hamiltonian, 
along with the expansion of the electric field in terms of 
(B6| and integrating over a time interval, [0,T], gives 
the amplitude for excitation after a time interval, T, 
C\,a,i.j,k.i(T). The rate of excitation of the detector, r ^i 
is given by the absolute square of this amplitude summed 
over all possible indices, and divided by the time inter- 
val, T. We take T to infinity, and obtain the result that 
would be expected from an application of Fermi's golden 
rule 



rw = 



2^2 



poo 

/ k, • Im \G{— k, x\\q, as || o , — w)] • k 

Jui/V 



dk 
2^ 
(25) 

Full details of the derivation of this result can be found 
in appendix [d| In general the transition rate (25 1 is non- 
zero, and it is therefore possible to extract energy from 
the centre of mass motion of a moving dielectric through 
coupling to the electromagnetic field within the dielectric, 
or in its vicinity |35j . This is a physical effect of the 



negative energy states within (20 1 that arises from our 
description of dissipation. 

Following the concluding discussion of section [II B[ it 
is clear that for moderate velocities we can only rely on 
(25) for low frequency excitations of the detector. For 



instance, if we were to consider a dielectric moving at 
lms^ 1 then an optical excitation oj ~ 10 15 Hz would cor- 
respond to integrating over wave-vectors, k > 10 15 m _1 
which is clearly far beyond a regime where the behaviour 
of the dielectric can be described with a susceptibility like 
9|. If the integrand in ( 25 ) decays rapidly (as it does in 
301)), then an uj in the GHz regime (k ~ 10 9 ) appears to 



be more appropriate. For relativistic velocities, (25) can 



be applied to much higher frequency excitations. Usually 
such velocities are completely inaccessible to experiment. 
However, they are worth considering in this case, for in 
experiments such as [5], a refractive index perturbation 
is created that moves through a dielectric at a fraction 
of the speed of light in vacuum. In the rest frame of the 
perturbation we have the dielectric moving at a uniform 
relativistic velocity that is subject to a weak stationary 
perturbation. It is therefore likely that the experimen- 
tally observed production of photons can be understood 



with some analogous formula to (25). This will be the 
subject of future work. 

The transition rate (25) vanishes when V = 0. This 



is because the argument of the delta function (D2) in 
the derivation of ( 25 ) is never zero when V 

r 



and so 

consistent with the concluding dis- 



>i = 0. However- 
cussion of section |II B| — the way it tends to zero depends 
on the response of the medium to large wave-vectors at 
a fixed frequency. If the atom is embedded within a di- 
electric medium this requires knowledge of the spatial 
dispersion of the susceptibilities. Yet when the atom is 



lm[K] 

V<c 



X 



V=c 



<i-Kr 



V>c 



Re[K] 



FIG. 2: As k is increased within the integrand of (25 I, the 



poles at ±Ko in ( |2Tp move along the red dotted lines in the 
directions indicated by red arrows. When V < c, Ko starts at 
a point on the imaginary axis indicated by a blue cross and 
moves vertically upwards. In the case when V — c, Ko starts 
at the origin. Therefore, for sub-luminal propagation, the 
poles never touch the real axis. Meanwhile, for super-luminal 
propagation, Ko starts at some point on the real axis, moves 
towards the origin, and then up the imaginary axis. In this 
case there is an imaginary contribution to the integral along 
positive Re[_RT] in (27) that comes from the need to negotiate 
the pole at Ko- 



hl a region of free space outside of a moving dielectric 
medium (see section IV 2 below), this is not necessary 
as the large wave-vector modes are very tightly bound 
to the surface. Therefore in this case the transition rate 
(f30l) falls to zero as V — > because the relevant modes 



do not reach the atom (see (30)). 



We now investigate ([25| in two simple cases; free space, 
and close to a surface. 



1. Free space 

Firstly, we show that To-vi vanishes in free space. The 
free space Green function is a function of the difference 
in the positions, a; & a;', and in k space is given by, 



G{k,u) 



"Mo: 



k ® k - (uj/c) 2 t 3 



(w/c) 2 (|fc| — uj/c — ir))(\k\ + uj/c + irf) 

(26) 

We transform this into the y-z plane of physical space, 
G(fc,JC|| - x||,w) = (2tt)- 2 J G(k, K ,u) exp[iK ■ (x\\ - 
x'^)]d 2 K, integrate over the polar angle in the K plane, 
and find the result, 



Im[G(fc,0,w)] =Mo 
KdK 



xlm 



— (-k 2 x (g) x 
J Q {K\x\\ -x'A 



V,|) + 1; 



2tt (K - Ko-ir]) (K + Ko + iv) 



(27) 



7 




electromagnetic field in the region of space outside of a 
moving dielectric. The equivalent case of atoms moving 
along surfaces has been investigated before (e.g. [TU1I5U] ). 
but usually the work is concerned with the drag force on 
the atom due to the vacuum field. The author is not 
being derived previously, or being 
states of quantum systems moving 



aware of equation ( 25 



applied to the interna 
over surfaces. 



For a detector situated outside of a moving dielectric 
(see figure [3]), we can obtain the expression for the Green 
function through taking ( 26 1 and writing it in the follow- 
ing form, 



FIG. 3: We consider the transition rate, To-n of a detector 
from the ground state, |0) to an excited state |1) as a function 
of the perpendicular distance, zo from a surface moving along 
the x axis at T = OK. The rate of excitation is given by ( 30 1, 



and decays with zo at some rate (grey dashed line) determined 
by an integration over the imaginary parts of the reflection 
coefficients of the moving surface. 



where Kq — \J (uj/c) 2 — k 2 , Jq is a zeroth order Bessel 
function, and Vy = d/dxu. Terms involving first order 
derivatives in x» have been set to zero given that they 
serve to replace the zeroth order Bessel function with a 
first order Bessel function, which vanishes as x\\ — > xi. 
The only possible contribution to the imaginary part of 
the integral in (27) would come from the need to ne- 
gotiate the poles if they meet the real line as 77 — > 0. 
However, when k > uj/c then Kq is imaginary (see figure 
[2]). In this case these poles never meet the real line, and 
the right hand side of (27) vanishes. Therefore, so long 



as V < c then (25) is zero in free space. This is consis- 



tent with existing findings that the vacuum is unstable 
to superluminal propagation (e.g. [29 ). 



2. Close to a surface 

We now consider the possibility of exciting an internal 
state of a system through interaction with the vacuum 



G(k, x» , x'u , Ld) 



d K i-nr-oii-a;,: 

(2n) 2 



|fc| 2 E A e A (fc)®e A (fc) 
(|fe| — w / c — irj)(\k\ + w/c + irj) 



(28) 



where the e\ t 2 represent the two unit vectors orthogonal 
to k, the unit wave-vector. We take a dielectric moving 
along x and lying in the x—y plane with free space in 
the region z > 0. The detector is positioned near to the 
surface (zq > 0). In this case when x' is set as the detec- 
tor position, only the first term within the integrand of 
(28) is non zero at the surface. In satisfying the bound- 



ary conditions at the interface we therefore ignore the 
term proportional to I3. Performing one of the K in- 
tegrals (say over k z ) enforces the dispersion relation in 
the remaining term, and ( 28 1 becomes an integral over 



incident propagating waves. The boundary conditions at 
z = can then be fulfilled with the usual reflection coef- 
ficients, rxy(k,u>), where A and A' each take the values 
1 or 2, with the meaning that polarization A is incident, 
and polarization A' is reflected (moving media mix po- 
larizations [31] ) 



dk y e *M«-»°) 

G{k,x h x\\ ,u)) = n \ I — 



ex(k, k v , isg) (g) e\(k, k y , is£,)e~ 



sf(2-Z0) 



^e\>(k,ky,i£) (g) e\(k,ky,-i€)r\x(k,ky,-i£,u})t 



-£0-Uo) 



AA' 



rM (2) (^|| -*||oU (29) 



where £ = Jk 2 + k 2 — uj 2 /c 2 and s = sign(z — z ). At x\\ = x\\ , the first term in the integrand of (29) becomes real 



due to quantities linear in s going to zero. Applying (25) to (29) then gives, 

,2 



rw = 



E 



dk f°° dk v -2g Z0 T 
2tt£ 



k ■ e\>(-k, k y ,i£,)n ■ e\(-k, k y , -i£)r\\>(-k, k y , -u>) 



ET ST fe-^| K .e A (-fc,^ ? 0| 2 Im 

y JUJ/V J-OO 2 ^ 



r\\(-k,ky,-i£,-w) 



(30) 
(31) 



where the second line is obtained through neglecting the 
mixed reflection coefficients, r\2 & r2i, which are a fac- 
tor of V/c smaller than r\x and r 2 2- Something like ex- 
pression (31) could have been anticipated from applying 



the reasoning of [6] which contains very similar expres- 
sions in the discussion of the phenomenon of quantum 
friction, based not on canonical quantization but on the 
fluctuation-dissipation theorem. For positive rest frame 
frequencies, w_ > the imaginary part of the reflection 
coefficient is positive. The range of the integration over 
k ensures that this is always the case. As the detector is 
moved away from the surface, the exponential decay on 
the right of (30) serves to reduce ro_n to zero. However, 



the dependence of this rate on the distance from the di- 
electric depends on the dispersion of the medium. This 
will be investigated in future work. 



V. CONCLUSIONS 

We have developed a formalism for the description of 
uniformly moving media that includes a full account of 



dispersion and dissipation. It was found that due to the 
presence of dissipation, the Hamiltonian describing the 
interaction between a uniformly moving dielectric and 
the electromagnetic field lacks a lower bound (alth ough 
this statement should be qualified as in section II B[ ) . Fur- 
ther to this it was shown that a stationary detector placed 
inside or outside a moving medium has a finite probabil- 
ity of excitation that depends on the specific dispersion of 
the medium, and in the case of being outside the medium, 
the distance from the surface. 
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Appendix A: Integral identity for susceptibilities 



To obtain expression (17) the following result must be applied, 
(a; 2 - n + n' + )[a EE (u) ■ a^.(w) + a EB (w) ■ < B M] 



dui 



(fi 2 _-rvy + )x BE (fi_) 



(w-Cl-- ir))(u) + fi_ + ir))(w - Cl' + - ir))(u) + + irj) (Q_ - - ir?)(H_ + Q' + + irj) 

(n' + 2 - n + n' + )x EE (n' + ) 



(fy + - - ir])(Q' + + + irj) 



(Al) 



along with similar expressions containing the susceptibilities, X EB , X BE , and % BB . The proof of (Al I is similar to that 
of the usual Kramers-Kronig relations. We first expand the left hand side of ( Al ) in terms of residues at poles and a 
principal part, 



du> 



{oj — f2_ — irj)(u> + + ir))(u) — Q', — irj){oj + Q,', + irj) 



(nl-n + n' + )im[x EE (n_) 



(fi' + 2 - n + n' + )im[x EE (n' + ) 



(f2_ - Q' + - ir)){Sl- + irf) (CL' + - fi_ - ir)){Sl' + + ft- + it}) 



duj 



- n + n> + )im[x EE (Lo)} 
(u 2 -n 2 _)(u) 2 -n> + 2 ) 



(A2) 



9 



where we have applied the result (see ffify) a EE • a^ E + a EB • a^ B = 2a;Im[% EE ]/7r. 

Now consider the following integral or a function, f(z), that is analytic within a closed contour C, 

1 f (z 2 -a)f(z)dz (z 2 1 -a)f(z 1 )-(z 2 -a)f(z 2 ) 



2ni J c (z - zi)(z - z 2 ) 



(zi - z 2 ) 



(A3) 



where a is real. If f{z) is analytic throughout the upper half complex plane, then C may be taken as the real line plus 
a semicircle at infinity, directed anticlockwise. If it is additionally supposed that f(z) is such that the integral over 
this infinite semicircle is zero (it goes to zero faster than 1/R, where R is the radius of the infinite semicircle), then 
( A3 ) may be written as, 



1 

2TTi 



(x 2 — a)f(z)dx 



(x — x\ — irj)(x — x-2 — iff) 



(x\ - op/pEi) - (x\ - a)f(x 2 ) 
(xi - x 2 ) 



(A4) 



where x\ & x 2 are real, and rj is as in the main text, an infinitesimal quantity that serves to shift the poles away from 
the real line and into the upper half plane. Expanding the left hand side of ( A4 1 into residues and a principal part 



1 



(x 2 - a)f(x)dx {x\ - a)f(xx) - (x\ - a)f(x 2 ) 



l-oo (X-Xi)(x-X 2 ) (Xl-X 2 ) 

we can obtain the following relation between the real and imaginary parts of f(z), 

x(x 2 — a)lm[f (x)]dx (x\ - o)Re[/(o;i)] - {x\ — a)Re[/(a; 2 )] 



(x 2 — x\){x 2 



4) 



(x 2 x 2 ) 



(A5) 



(A6) 



when (A6| is applied to (A2) we obtain (All. 



Appendix B: Field commutation relations 



Here we show that the substitution given in ( 18 ), plus the assumed commutation relations between C & ( 19 1 are 



consistent with the commutation relations that must hold between the canonical momenta and the field amp 
In the gauge where V • A = 0, these canonical commutation relations are given by [2"2"] . 



itudes. 



A(x,t),tl A (x',t) 
where S±(x — x') is the transverse delta function, 

d 3 k 



ihd±(x — x') 

iht 3 S(uj - Lu')6 i3) (x - x') 

ih± s 5(u - lu')S (3) (x - x') 



(Bl) 



8 ± (x -x')= f d k - (k 2 t 3 -k®k) e ik < x - x ') 
V ' J (2^) 3 fe 2 V 7 



(27r) 3 fc 2 

The vector potential operator is expanded in terms of the C\ operators as follows, 



A(x,t) = Y, 



dk 
27 



d 2 x[ 



h r 



f\(k, x\\ , x\ , n) • C x {k, x'« , 7 ft)e 4(fex - (n+yfe)t) + h. 



(B2) 



with similar expressions holding for the remaining field operators in (Bl|. To find these expansion coefficients, one 
must solve the classical problem of section II A and then perform the substitution (18). Inserting the operators in the 
form given by (B2) into (Bl) we obtain the conditions that the expansion coefficients must satisfy in order that the 

ik( - x - x,) f x A (k, x|, , s^fi) • (fc, x\ { , xm,n) = iS x (x - x') 



canonical commutation relations be satisfied 

dVL 



E 



dk 

2^ 



E 



dk 
2^ 



d SB 111 



d 2 x\\ 1 
dQ 



(B3) 



ik(x-x') 



o 2 7 fi 



fx M (k,x\\,x\\i,ty ■ /nt (k,x[,x n ,ri) - c.c. 



it 3 S(u}-u}')S^(x-x') 

(B4) 
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where we have applied, f A (—k,x\\,x'n,—Cl) = f A (k,x\\,xi,£l) and a similar relation for fn A (k, xy,Xy, fi). In this 
case the expansion coefficients for the electromagnetic field operators are given by, 



f A (k, xy fi) = G±(k, £C|| , x n , fi + Vk) ■ O x (k, 

/|(fe, xii ,X|h, fi) = i(fi + Vk)G(k, xii , xp.fi + Ffc) • 6 A (fe, xihjfi) 



The electric field expansion coefficients (B6| are here given for reference, and in the above we define 

d A (fc,cC||,fi) = [i(n + Vk)al x (x hl n) + a^(x||, 7 fi) • Vx] 
6 A (fc,xy,fi) = [-i(n + Vk)a EX (xn,yQ) + x V ■a BA (x,,,7fi)] 



(B5) 
(B6) 



(B7) 



where, as in the main text we have V = ikx + Vy and the operation, Gx V is equivalent to the curl with respect 
to the right hand index: Gx V= e^c^-G^. The Green function within the above formulae is a solution of 



V X 



/x _1 (x|| , fi_) • V X GO, xy , x\ , fi) - fi 2 e(cc|| , fi_) • G(fc, xy , x\ , fi) 

+ ifi[x BB (xi|,fi_) • V x G(fc,xy, x! h fi) - V x (x BE (xy,fi_) ■ G{k,x h x\,Q))] = l 3 <5 (2) (xi| - x[) (B8) 



The susceptibilities are defined as in (|9|), with e(xy,fi_) = e l3 + X EE (x|| , fi-), and fi 1 (x\\,Q-) = /i l ta — 
Xbb( x \\ > fi-)- We note that in a moving medium the Green function does not satisfy the reciprocity condition, 



G(x,x',fi) ^ G 1 (x',x,fi) 



(B9) 



and therefore G does not satisfy ( |B8[ ) with respect to its second argument, as is usually the case. Instead it satisfies 
(|B8| with V -> V ( Xeb -> -X«bM) 

G(fc,xj |) X||,fi) X ^* • 1 (ecu , 0_)1 x V* - fi 2 G(A;,x||,X||,fi) • e(x||,fi_) 

-iSl[G(k,xf l{ ,x\\,Q) X ^•x BE (x| | ,fi-)-(G(fc,x; | ,X||,fi).x EB (^||>^-)) X V*] = l 3 <5 (2) (xy - x\ { ) (BIO) 
The transverse Green function which appears in the expansion of the vector potential operator is related to ( |B8[ ) by 

Gj_(k, xi|,XM,fi) = / d 2 xi|i<5j_(fc,X|i — X|h) • G(fc,X||i,Xn,fi) 



( BU ) 

with, Sj_(k, X|| — xj|) = J" <fa;<$i(x — x')e~ tktyX ~ x \ The following result will be important in proving that the commu- 
tators take the correct form, 

/ d 2 x n G(k, xy ,x n ,n + Vk)-J2 o[(k, xy!, fi) • 6 A (A, xy!, fi) • G f (fc, x\ .acp.n + VA) 

^ A 

r G(/fc,xii,Xn,fi + yfc) -G f (fc,a;n,a;||,fi + Ffc)j (B12) 



This can be proved from taking (BIO) multiplied on the right by G\ subtracting the adjoint of (BIO) multiplied on 
the left by G, and integrating over xy. This can then be seen to equal (B12) after and application of (B7) and (p" 



Throughout these manipulations we must take care to remember that the symbol V contains a component ikx that 
reverses sign relative to Vy when integrating by parts. Explicitly if we take two vector fields, V(k, xy) & W(k, xy) 

then, neglecting boundary terms, J d 2 x\\V ■ V X W = J d 2 x\\W • V* X V = j d 2 x\\V X V* • W. 

The remaining expansion coefficients are as follows: the canonical electromagnetic field momentum expansion 
coefficient is given by 



(fc,xy,X|h,fi) = -<5 (2) (xy - X|h)a EX (x|h,7fi) 



i(fi + Vk)e(x {] , 7 fi) + X*b(*\\ , 7^) • V X 



G(fc, xy , X|h , fi + Vk) • O x (k, X|h,fi) (B13) 
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the expansion coefficients for the oscillator field operators are 

,x\ {1 ,n) = M AE <5( 7 ^ - «)*( 2 >(x|| - xp) 

[i(fi + 7fe)aJ E (x||,w) + a^ E (x||,w) • Vx] • G(fc,aj||,aj|| 1 ,0 + Vk) ■ &\{k,x\\i,Si) 

{ui — 7fi — irj)(u) + 7f2 + £77) 

and the expansion coefficients for canonical momentum operator for the oscillator field are 

(jfe,aj||,xi|i,n) = -i7 2 n/x (fe,a:||,a5||i,fi)- 



(B14) 



(B15) 



We now insert | [B5| , dB13] >, |B14| & pl5) into p3|]B4| . In the case of the canonical commutation relation for the 
electromagnetic field (B3) this gives, 



dk 
2n ( 



ik(x- 



2tt 



G±(k,x [l ,x' n ,£l + Vk)- 



-i(fl + Vk)e(x[ , Tfi) + X V . Xbe (x(| , Tfi) 



-i(n + VA)e*(4,7fi)+x V . % * E (x' 7 ft) 



= 5l(x-x') (B16) 



where we have applied ([9]) and (B12). The adjoint of the Green function, G^ in (B16) is transverse with respect to 
the left hand index. 

The Green functions and the susceptibilities are both analytic functions of frequency in the upper half plane, with 
the complex conjugates of these quantities being analytic in the lower half plane. Therefore we can split the integrand 
of (B16) into two pieces — one analytic in the upper half plane, and one in the lower half. We can then deform the 
integration contour over frequency from being along the real line, to following a semicircular path from — oo to oo, 
in either the clockwise (C + ) or anticlockwise (C~) sense, depending on where the integrand is analytic. Along C we 
have, |f2| — ¥ oo, where the Green function can be taken to satisfy the free space equation and is therefore equal to 



G(fe,X||,x(|,Q) -> -/x 



d 2 K 



^i 3 -fc<g,fc \ 

« 2 fSV /„2\ 



(B17) 



where h = kx + K. As fc is a real vector, the deno minator in ( |B17[ ) is non-zero over the entire contour. Therefore, 
G(fc,X||,Xn,ft) ->• -^^l 3 <5 (2) (a;|| - xf,), and (B16l simplifies to, 



dk 
2? 



r dn 


r dn ' 


? 


Jc+ 4 


J c - 2vrSl 





iS±(x — x') 



(B18) 



writing Q, — |f2|e je , and performing the two contour integrals in the square brackets gives a factor of i, and the equality 
is fulfilled. Therefore the expansion of the electromagnetic field operators in the form (B2) with the coefficients (B5) 
and (B13) is consistent with the canonical commutation relations. 

For the oscillator fields, inserting ( |B14 ) and ( |B15 ) into (B4|, we find that after applying (B12), all terms cancel in 
the product [/ ^(fc, x\\, xyx, 0) • / n ^. ^ (fc, xi, xyx, 0) — c.c] except for the first involving products of delta functions. 
We are then left with 



E 



dk 
2^ 



d xih 



dn r 
27ft 



rfne ik ( x - x 'h 3 5 XE 5 XE 6(~fn - u)S(jQ - uj')5 {2) {x« - xyx)£ (2) (x|, - xim) 



x - x) (B19) 



Performing the integrals and summation on the left of ( B19 1 shows that the equality is fulfilled, and this completes 
our demonstration that ( 18 1 is a canonical transformation. 



Appendix C: Including the centre of mass as a dynamical variable 



A natural reaction to a Hamiltonian for a quantum field that includes negative energy excitations, such as (20 1, 
is that it is physically unacceptable. If the Hamiltonian included all of the possible dynamical degrees of freedom 
as canonical variables then this would probably be true. However, in our case the centre of mass velocity of the 
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dielectric medium is assumed to be fixed. We might therefore interpret the negative energy as an accounting device 
for the external energy input required to keep the motion constant. To show that this interpretation is reasonable, 
we now demonstrate that including the centre of mass as a dynamical variable removes the terms that cause ( 20 ) to 
be unbounded from below. 

We start from ([!]), adding in an extra term to represent the kinetic energy of the centre of mass. The centre of 
mass is only a meaningful concept to first order in V/c, and we therefore assume a motion of the medium such that 
we can set 7 ~ 1. For arbitrary motion we should use the centre of energy, and the treatment becomes much less 
straightforward. To this order the total Lagrangian is, 



L = -MV 
2 



d X [^2p H - ^t - ^int] 



For these purposes, we assume that V = Vx as in the main text, and construct the Hamiltonian as in section II B 
The canonical momentum associated with the centre of mass is, 



dL 
W 



/poo 
d 3 x J duj [V (g) X u ■ U Xul - a EB (x\\,u)X u X B + a BB {x\\,u)Y u X E/c 2 ] 



and therefore the Hamiltonian is, 



H = -MV 1 



(E 2 



2 B 2 ) + 1 jf° doj [n^ + + lo 2 (Xl + Y 2 U ) - a BB (u})Y u .(B + VX E/c 2 )] J 

(CI) 

The term, V ■ (V <S) X u ) • TLx„ present in (15) has been removed from (CI) through the inclusion of the centre of 



mass dynamics. As explained in the main text, it is this term that is responsible for the lack of a lower bound to the 



Hamiltonian. It is not possible for (CI) to be decreased to an arbitrary negative value. 



Appendix D: The transition rate of a detector interacting with a moving medium 



To find the transition rate, we first find the probability amplitude associated with the atom making a transition 
into an excited state, given that the field is also in some given state. In (24) it was found that the rate of change of 
this amplitude is given by 

(\,a,l, m ,n,p{t) = ~(l\,<r,l, m ,n,p\K ' E(x ,t)e iult \0) 

where the matrix element associated with the atomic transition has been calculated. To calculate the matrix element 
associated with the combined system of th e electromagnetic field and the dielectric medium we use the expansion of 
the electric field in the form given in |B2| . We then expand the C\ operators in terms of the c a \ operators using 



The orthonormality of the number states given in ( 23 ) is then applied to give 

d 2 xi 



1 



dk 
2^ 



dn 



i[-kx + (n+Vk+Lo)i\ 



Integrating this over a time interval [0, T], squaring the resulting expression and summing over all possible final states 
then gives 



J2 \Cx,a,l, m ,n, p (T)\ 2 = — d 2 x\ / *||o,*j|,n) 
i Bm „ J J J \ 



\,a,l,m,n,p 



2 2sm 2 [(n + Vk + uj)T/2] 



(Dl) 



where we assume that the value of the expansion coefficient is equal to zero at t = 0, and we have applied the 
completeness relation (22 ). Dividing (Dl ) by the time interval, T gives the transition rate, ro_n. As T — > oo, we can 
apply the identity 



Asin 2 [(n + Vk + uj)T/2} 
t^o {n + Vk + uj) 2 T 



2n8[{Vl + Vk + u))] 



(D2) 
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and the transition rate simplifies to 

rw = - 



l r-u/V r 2 

— / dk d 2 x {l Y^ — k- f*(k,x ll0 ,x {l ,-u-Vk) 



(D3) 



Inserting (B6) and applying (B12) then gives the transition rate in terms of the Green function, 

2uj 2 



rw = -- 



iv 



dk G(k,x\\ ,x\\ ,u)) - G'(k,xu ,xu ,uj) 

— — K • " - " " • K 

2tt 2% 



(D4) 



When a medium is reciprocal then it satisfies (B9), and G T (k, x»q, xuq, lo) = G(k, x\\q, snyo, w). As stated in appendix 
fBl for moving media (B9| is not satisfied, and we might expect this is affect ro_>i. Yet apparently the non-reciprocity 



has no effect on the transition rate of interest here (|D4j), for we contract the indices of the Green tensor with the 
symmetric tensor k ® k, which gives 



r, 



2co 2 



0-s-l 



dk 



2u? 



iv 



k ■ Im[G(fc,a3|| ,a;||o,w)] • k 



iv 



dk 
2^ 



K • Im[G(— k, xyo, aiiio, — u>)] • k (D5) 



which is the expression given in the text (25). 
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